Enhanced motility of a microswimmer in an elastic tube 
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We analyse the effect of confining elastic boundaries on the motility of a model dipolar microswim- 
mer. The flexible boundaries are deformed by the velocity field of the swimmer in such a way that 
the motility of both extensile and contractile swimmers is enhanced. The magnitude of the increase 
in swimming velocity is controlled by the ratio of the swimmer and elastic timescales, and the dipole 
moment of the swimmer. We explain our results by considering swimming between inclined rigid 
boundaries. 



Introduction:- Confinement in elastic and rigid envi- 
ronments is a key concept affecting fluid transport and 
locomotion in microscopic systems, ranging from molec- 
ular motors to single and multicellular self-propelled or- 
ganisms The low- Reynolds number world of micro- 
swimmers is often crowded by passive and active bound- 
aries that can be exploited by swimmers to enhance their 
motility, or as defence mechanisms which act to hinder 
their activity. 

Not surprisingly, the locomotion of pathogens close to 
flexible and rigid surfaces can play a key role in their 
success or failure as infectious agents. Motile pathogenic 
bacteria and parasites move through the extracellular 
matrix of the host to invade the vascular system [2[ and, 
in the case of neurodegenerative infections, even cross 
the blood-brain barrier |3[. For example, pathogenic 
sphirochaetes responsible for syphilis Q , leptospirosis Q 
and Lyme disease [6j, have been reported to swim in 
microvasculature channels Q, renal tubules [H[, and to 
invade the intracellular junctions between endothelial 
cells • Defensive mechanisms against infections also in- 
volve close interaction between complex elastic surfaces 
and pathogens. Flexible-surfaced leukocytes prey on 
motile microorganisms [8] , such as E. coli and P. aerugi- 
nosa, and the viscoelastic gel-like lining the stomach acts 
as a barrier against H. pylori [9[ , a flagellated rod-shaped 
bacterium linked to chronic gastritis and stomach cancer. 

Understanding the way in which microswimmers be- 
have in confinement can also lead to novel applications in 
microfluidics and biotechnology. Recent advances in flow 
imaging and cell tracking have made it possible to mea- 
sure the scattering function of sperm cells and algae from 
surfaces, encouraging clever rectification applications in 
microfluidic chambers 

[3 EH, 

and the increasingly ac- 
curate magnetic steering of artificial micro-robots inside 
microchannels [12|, [l3j has promising applications in non- 
invasive drug-delivery technologies. 

Motivated by these observations, in this Letter we for- 
mulate the following questions: what is the effect of the 
activity of a microscopic swimmer on a bounding elastic 
surface and, what is the effect of the interaction with the 
surface on the motility of the swimmer? While previous 
efforts have mainly focused on swimming strategies that 



exploit rigid and elastic surfaces to overcome the scal- 
lop theorem f3, EH , the effect of confinement by rigid 
and flexible boundaries on the motility of microswim- 



mers remains largely unexplored [16|, LL7j . By providing 
a microscopic footing, single-particle interactions can be 
used to build up coherent long- wavelength hydrodynamic 
models of active matter subject to viscoelastic hetero- 
geneity [13, LL9| and confinement [2CH22| to model, e.g., 
peristaltic pumping and cytoplasmic streaming [23|, |24| . 

Our results are summarised in Fig.Q] which shows sim- 
ulations of a model swimmer moving along the centre line 
of an elastic channel, following the direction indicated by 
the thick arrows. The swimmer, which has advection 
timescale r s , is force- free and acts as a force dipole on 
the fluid, of strength p. The resulting flow pattern, pro- 
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FIG. 1. (Colour online) (a) and (b) Simulations of dipolar 
extensile and contractile swimmers passing through an elastic 
tube. The deformation of the channel, protruding on average 
towards extensile swimmers and away from contractile ones, 
is set by the competition between the swimmer velocity field 
(projected planes) and the elastic response of the boundary, 
(c) and (d) The increased friction coefficient due to the solid 
boundaries leads to a net increase in velocity, Vf, which is 
larger for extensile swimmers. The amplitude of the response 
for swimmers with the same self-swimming speed v s depends 
on the ratio between swimmer and elastic timescales, r s /rf, 
and the swimmer dipolar strength, p (So ~ p). 
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jected in the bottom of each frame, can be extensile (fluid 
is pushed out from the ends of the swimmer and drawn 
in to the sides) as shown in panel (a) or contractile (fluid 
is pushed out from the sides and pulled in to the ends) as 
shown in frame (b). The shape of the flexible boundaries 
is set by a competition between the activity of the swim- 
mer, which tends to deform the tube and is characterised 
by r s and p, and the elastic response of the confining 
boundaries, which resists such deformations and is con- 
trolled by the bending relaxation timescale rp . Note the 
different shape deformations for extensile and contractile 
swimmers. 

The effect on swimmer motility is presented in frames 
(c) and (d), which show the change in swimmer speed, 
relative to the free swimming speed, resulting from the 
interaction between the swimmer and the boundary. Our 
simulations show that both kinds of swimmers move 
faster through the channel, but extensile swimmers ex- 
perience a stronger speed-up than contractile swimmers. 
Moreover, the variation in motility is strongly dependent 
on the ratio of swimming and elastic time scales, r s /rf, 
and on the dipole strength. 

To explain these results we argue that the interaction 
between the model dipolar swimmer and the boundaries 
can be understood in terms of shift and tilt deformations 
to the boundary. We show that shift deformations (to- 
wards the swimmer) always increase the swimmer speed 
by an amount that is proportional to the self-swimming 
speed and independent of the strength or direction of 
the velocity field created by the swimmer. Conversely, 
tilt deformations, corresponding to a local inclination of 
the boundary, couple to the particular swimming stroke 
and can result in enhanced or reduced motility depending 
on the swimming pattern. Using these results we are in 
a position to argue that the hydrodynamic coupling be- 
tween the swimmer and a flexible boundary always leads 
to channel deformations which favour the passage of the 
swimmer. 

Model swimmer:- Locomotion at low Reynolds num- 
ber relies on swimming strokes that are non-reciprocal 
in time. Many microorganisms achieve this by deform- 
ing their shape in such a manner that the local drag 
coefficient varies along their bodies. We coarse grain 
such a feature by considering a model swimmer made 
up of anterior and posterior spheres, a and p respec- 
tively, joined by a link of prescribed but variable length 
/ [25| (see Fig. 15(a)). The local drag on each sphere 
follows from the Stokes Law, = <^v^, where Vi is 
the velocity of the z-th sphere and & its friction coef- 
ficient. The swimmer is subject to the force- free con- 
dition, F p + F a = 0, and to the kinematic constraint 
dictated by the swimming stroke, v a — v p — Z, where the 
dot implies differentiation with respect to time, t. As 
a consequence, the instantaneous speed of the swimmer 
obeys v = (v p + v a )/2 = Z(f p - f a )/2(f p + f a ). 

For free, unbounded, swimmers locomotion can only 



be achieved if the friction coefficients vary in time, intro- 
ducing a non-reciprocal deformation of the body of the 
swimmer. This can be readily verified by considering the 
Stokes drag coefficient & = 6^77^, where 77 is the vis- 
cosity of the fluid and ai is the instantaneous radius of 
the z-th sphere. Without loss of generality, we can set 
I = lo + 51 sin(co>£), a p = ao + Sa sm(ojt + Aip + Sip), and 
&a = ao+5a sm(ujt+Aip — 5ip) to obtain the average speed 
of the swimmer over one stroke, v s = f^_^° dtv/(27r/uj) = 
SlojSasmSip cos Aip/4a + O(5a/a ) 2 |26|. 

The instantaneous forces exerted by the swimmer on 
the fluid can be averaged to obtain the net flow induced 
by the swimmer. Due to the force- free condition the far- 
field velocity field is dipolar, (u) = p(3 cos 2 <p— l)r/87T77r 2 , 
where r = rr is the displacement vector from the centre 
of mass of the swimmer to a point in the fluid, and <p 
measures the angle from the axis of the swimmer to r 
in the counter-clockwise direction. The dipole strength 
p = 37TT]5luj5alo cos Sip sin Aip/ 2 sets the typical magni- 
tude of the field, and its sign, controlled by the phase 
differences, Aip and Sip, determines whether the flow pat- 
tern is extensile (p > 0) or contractile (p < 0), depicted 
in frames (c) and (d) of Fig. [TJ 

Rigid boundaries:- We first focus on the interaction 
between the model swimmer and a bounding solid sur- 
face. We consider the situation where the swimmer and 
the surface are aligned along a single axis of revolu- 
tion, as shown in Fig. [2(a). For such configurations 
the drag coefficient increases, and can be written as 
& = 67rr]ai/(l + C[ a i/hi])i where ((x) < is a correction 
term accounting for the additional friction offered by the 
solid. In general, £ is a decreasing function of ai/hi, and 
its functional dependence is specific to the geometry of 
the surface. The interaction with the solid will lead to a 
net contribution to the swimming speed, v w . For planar 
walls (hi = ho) and prescribed deformations, the swim- 
mer 'grips' the surface, exploiting the higher resistance 
offered by the walls. More specifically, 

- = - »-(£)nt (1) 

where (0 = C[ a o/^o], etc. The scaling of Eq. (pQ) indicates 
that the swimmer always moves faster in confinement 
(Co < 0), by a factor controlled by the proximity to the 
wall and its hydrodynamic resistance. A similar speed- 
up effect has been obtained for a waving sheet swimming 
between flat walls [16| . 

Note that Eq. (pQ) predicts that both extensile and 
contractile swimmers increase their speed by the same 
amount in flat confining geometries. This is because the 
hi do not vary along the swimmer trajectory. This, how- 
ever, changes if the walls are inclined relative to the swim- 
mer path. As illustration, consider small wall inclinations 
such that hi = ho + Si. This will lead to an additional 
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FIG. 2. (Colour online) Swimming in confinement, (a) Dipo- 
lar swimmer composed of two spheres, a and p, in the middle 
of planar, inclined walls, (b) Extensile swimmers push on 
the posterior end of their bodies (black arrows), and experi- 
ence a speed-up when in diverging channels, light red arrow, 
(c) Conversely, contractile swimmers pull on their anterior 
end during the power step of their stroke (black arrows) , and 
their speed is increased in converging channels (light green 
arrow), (d) and (e): Contour plots of the normalised contri- 
bution to the speed of extensile and contractile swimmers due 
to inclined boundaries, Vi/v s , as a function of the average dis- 
placement, Ai, and net amplitude, A2, of the channel. The 
dotted paths correspond to the variation of Ai and A 2 along 
the top curves in Figs. [TJc) and[TJd). 

contribution to the speed of the swimmer, 

»—te)( a - +a ^) z « +o<a ^- < 2 » 

The first term, which reflects the net displacement of the 
position of the surface from the reference value, ho, is 
controlled by the shift, Ai = (S p + 5 a )/2ho. The sec- 
ond term results from the local inclination of the surface, 
measured by the tilt, A 2 = (S p — S a )/2ho. Both terms 
scale with the function Z = Z[ao/ho] < 0, which en- 
codes the strength of the resistance provided by the wall 
through Z = [C - (a /h ) (C' 2 /(i + C) - C")] /(I + 0- 

While the first term in Eq. (j2]) can be absorbed as a 
simple offset of ho in Eq. (pQ), the A 2 -term introduces an 
interplay between the local structure of the wall and de- 
tails of the swimming stroke, as indicated by the phase 
dependence. Both extensile and contractile swimmers 
achieve propulsion by pushing on their posterior end dur- 
ing the expansion of the link (a p > a a ), and by pulling 
on their anterior end during contraction (a a > a p ). 
For extensile swimmers, as depicted in Fig. Efb), the 
'power step' corresponds to the extension step, where 
the spheres have relatively large radii. Diverging con- 
strictions, A 2 < 0, amplify the strength of this step, 
leading to an additional speed-up. For converging con- 



strictions, however, the power step is weakened by an 
amount controlled by A 2 , which can lead to a net slow- 
down. Similarly, contractile swimmers generally swim 
faster in converging constrictions-as shown in Fig. El^c)- 
and slower in diverging ones, relative to the planar wall 
reference configuration. 

The contour plots in Figs.[2fd) and[2fe), show how the 
interplay between wall proximity and wall inclination can 
lead to both speed-ups and slow-downs for the swimmer. 
Note that v\ can be of the same order of magnitude as 
v s . To obtain these plots we have set Zo according to 
Faxen's correction to Stokes Law [27], whereby ((x) = 
Ax + Bx 3 + Cx b , with A = -1.004, B = 0.418 and 
C = —0.169. This particular friction law is valid for a 
pair parallel plates, relevant for swimmers confined in 
microfluidic chambers, and is expected to hold for our 
treatment in the limit of gently, locally, inclined surfaces. 

Elastic boundaries:- We are now in a position to ex- 
plain the results in Fig. 1, showing that both extensile 
and contractile swimmers move more quickly in elastic 
channels for a prescribed swimming stroke. However 
first we summarise details of the simulations that led 
to these results. The swimmer moves along the cen- 
tre line of a permeable cylindrical boundary formed of 
eight long bead-spring filaments. Each of the filaments is 
composed of forty spherical beads of radius <2f joined by 
Hookean springs of potential U s = k s (rij — r eq ) 2 /2, where 
Tij = Tj —Ti is the displacement vector between the i— th 
and j— th beads, r eq is the rest-length of the link and 
k s its elastic constant. We consider inextensible links, 
which we accomplish by setting the stretching relaxation 
time, T s tretch ~ ??af//c s , to be the fastest timescale in the 
system. Triplets of sequential beads, labelled ijk, are 
coupled by a bending potential Uh = —k^r^ • Yjk/rijVjk, 
with bending constant fcb- Both potentials are consid- 
ered in calculating the total force acting on each bead, 
Fi = — V(L / "b + U 8 )i. We consider the hydrodynamic in- 
teractions arising from the forces exerted by all beads 
on the fluid, whereby the velocity of a particular bead 
obeys vf = Hf? F? , where Greek indices denote Carte- 
sian components and repeated indices imply summa- 
tion. The mobility matrix is considered up to first or- 
der in aij corresponding to Oseen's approximation, 

H?f = SapSij /6tt W + (1 - SijftSap + rg-rg/r?.]/87rrp-y, 
where Sij is the Kronecker delta. Simulations are run 
by considering the passage of a single swimmer through 
initially fiat filaments and integrating the equations of 
motion using Euler's algorithm. 

The enhanced motility observed in Fig. [T] can be ra- 
tionalised in terms of the local deformation of the chan- 
nel. Qualitatively, it is apparent from Fig. [3] that ex- 
tensile swimmers deform the boundary to create a lo- 
cally diverging channel, while contractile swimmers in- 
duce a local converging geometry. These configurations 
lead to the expected speed-up, according to our argu- 
ments. Two main effects set the deformation of the fila- 
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FIG. 3. (Colour online) Typical deformations of an elastic 
boundary caused by (a) extensile and (b) contractile swim- 
mers. Due to the flow-field pattern, shown as light coloured 
arrows, the interface is deformed to a locally diverging or con- 
verging shape for extensile and contractile swimmers, respec- 
tively. The local deformation of the filaments is quantified by 
interpolating S p and 5& at the x— coordinate of the respective 
beads. 



merits, and in turn the back flow Vf. On the one hand, 
the filament is deformed by the velocity field of the pass- 
ing swimmer, of typical magnitude ~ p/v^o^ over tne 
swimmer self-propagation timescale r s = lo/v s . On the 
other, the elastic response of the filament resists defor- 
mations, and can be quantified by the bending relaxation 
timescale rp ~ 7?afA 3 /G, where A ~ ho is the wavelength 
of the perturbation and G ~ k^r^ is the bending mod- 
ulus. The growth-rate of perturbations to the bound- 
ary thus obeys, S « vo — S/r^ which suggests a scaling 
5 ~ S [1 - exp(-r s /r f )] /(r s /r f ), where S = vor s is the 
typical amplitude for freely-deformable chains. To ex- 
plore the interplay between these effects in detail, we 
ran simulations considering a pair of initially parallel 
filaments (Fig. [3]) and explored a wide range in r s /rf. 
(This simple configuration only reduces the magnitude of 
the back-flow with respect to the many-filament simula- 
tions, and does not alter results qualitatively.) Figs.QJc) 
and [1(d) show the dependence of Vf on r s /rf for exten- 
sile and contractile swimmers of identical free-swimming 
speed v s but with different dipole strength, p. As sug- 
gested by the simple scaling argument, curves tend to 
reach the rigid-boundary limit at r s /rf ~ 1, when the 
rigidity of the filaments suppresses deformations. Ad- 
ditionally, the amplitude of the curves scales with 5o, as 
expected for swimmers with stronger dipole moment that 
can induce larger deformations to the boundary. 

The particular shape of the curves shown in Fig. [T] 
is set by the subtle interplay between the local inclina- 
tion and shift induced by the swimmer on the elastic 
boundary. To illustrate this point, we measured the lo- 
cal amplitudes 5 P and S a as indicated in Fig. [3] for the 
top curve in Figs. [1(c) and H(d). We then calculated 
the corresponding Ai and A2 values, which we depict 
as paths in (Ai, A2)— space in Figs. [2(d) and [2(e). Ex- 
tensile swimmers tend to pull on the boundary, shifting 
it to closer positions and to higher inclinations simulta- 
neously. Both effects decrease with increasing filament 
rigidity, leading to a trajectory in (Ai, A2)— space that 



runs monotonically from negative values of both parame- 
ters to the origin, and is reflected in the smooth decrease 
of V{ with r s /rf shown in Fig. [1(c). Contractile swimmers, 
on the contrary, strongly deform the boundary, creating 
high local inclinations, but only shift its position weakly. 
For intermediate rigidity of the filament, while the in- 
clinations die out, the boundary is shifted to a closer 
position to the swimmer. This contributes to speeding it 
up. However, for higher rigidity the shift decays and this 
effect vanishes. This competition leads to a turning tra- 
jectory in (Ai, A2)-space which shows a good correlation 
with the non-monotonic behaviour observed for Vf —the 
maximum in the top curve of Fig. [1(c) corresponds to the 
right-most point of the path in (Ai, A2)-space (red dot) 
in Fig. [2(e). 

We have described the hydrodynamic coupling between 
the self-propulsion of microscopic swimmers and the re- 
sistance offered by elastic confining surfaces. Both con- 
tractile and extensile swimmers deform the surfaces in 
such a way that they are speeded up. A similar mi- 
croscopic mechanism is likely to underlie the enhanced 
motility observed in viscoelastic fluids [l|. We have also 
shown that there is an asymmetry in motility in locally 
inclined constrictions with extensile (contractile) swim- 
mers moving faster in diverging (converging) channels. 
This asymmetry has potential as a means of rectification 
of microorganisms according to their swimming stroke, or 
to create hydrodynamic barriers that prevent the entry 
of pathogens into permeable liquid reservoirs. We hope 
that our research will motivate experiments to explore 
these possibilities. 
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